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Abstract. The study by White (1984) on the growth of angular momentum in 
dark haloes is extended towards a more detailed investigation of the spin parameter 
A = Ly/E/GM 2 - 5 . Starting from the Zel'dovich approximation to structure formation, 
a dark halo is approximated by a homogeneous ellipsoid with the inertial tensor of the 
(highly irregular) Lagrangian region T from which the dark halo forms. Within this 
approximation, an expression for the spin parameter can be derived, which depends 
on the geometry of T, the cosmological density parameter i? , the overdensity of the 
dark halo, and the tidal torque exerted on it. For Gaussian random fields, this expres- 
sion can be evaluated statistically. As a result, we derive a probability distribution of 
the spin parameter which gives A ~ O.OTl^Qg, consistent with numerical investigations. 
This probability distribution steeply rises with increasing spin parameter, reaching its 
maximum at A ~ 0.025. The 10 (50,90) percentile values are A = 0.02 (0.05,0.11, re- 
spectively). There is a weak anticorrelation of the spin parameter with the peak height 
v of the density fluctuation field A oc v~ - 29 . The dependence on Qq and the variance a 
of the density-contrast field is very weak; there is only a marginal tendency for the spin 
parameter to be slightly larger for late-forming objects in an open universe. Due to the 
weak dependence on a, our results should be quite generally applicable and independent 
on the special kind of the fluctuation spectrum. 
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1 Introduction 

We adopt the widely accepted view that galaxies acquire angular momentum from the 
tidal torques of surrounding matter (Hoyle 1949, Peebles 1969, Doroshkevich 1970, White 
1984). Most of the angular momentum L is gained during the linear phase of the collapse, 
while an overdense region is expanding slower than the Hubble flow (i.e., collapsing in 
comoving coordinates) and angular momentum is growing linearly with time (L oc t, 
White 1984). After the turnaround point is reached and the collapse becomes nonlinear, 
further gain of angular momentum is strongly suppressed because it becomes increasingly 
difficult for the tidal field to act on the shrinking radius of the forming dark halo. Nu- 
merical simulations by White (1984) support this picture for the growth of the galactic 
angular momentum at least in a statistical sense. Numerical simulations by Efstathiou 
& Jones (1979), Barnes & Efstathiou (1987) and Warren et al. (1992) showed an average 
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spin parameter A = Ly/E/GM 2 - 5 of about 0.05... 0.07, where L, E, and M are the 
total angular momentum, energy, and mass of the dark halo, respectively, and G is the 
gravitational constant. Only 10% of the dark matter haloes in such simulations have a 
spin parameter of either less than 0.02 or more than 0.11. Such low values for the spin 
parameter indicate that the dark haloes are supported by the velocity dispersion of their 
constituents, and that the flattening of the haloes is a consequence of an anisotropic ve- 
locity dispersion and not of rotation. Furthermore, the spin parameter seems to be nearly 
independent on the peak height and the spectral index of the fluctuation spectrum. 

Up to now, several analytical studies investigated the growth of angular momentum 
in hierarchical scenarios, most of them within the CDM cosmogony in i?n = 1 model 
universes (Hoffman 1986, 1988, Ryden 1988, Heavens & Peacock 1988, Quinn & Binney 
1992). These investigations yield spin parameters on the order of 0.07, consistent with 
iV-body simulations. The more detailed studies of Ryden (1988) and Heavens & Peacock 
(1988) employ the peak formalism of Bardeen et al. (1986). Furthermore, Hoffman 
(1986) argued that the average spin parameter should be anticorrelated with the peak 
height of density fluctuations. Based on dimensional arguments, White (1994) suggested 
that the spin parameter in an open universe (!2o < 1) should be systematically lowered 
by f2 o e (t), where fi(t) has to be taken at turnaround time. However, a more systematic 
study of the spin parameter and its distribution remains to be done. 

In this paper, we investigate the statistical properties of the spin parameter in terms 
of global properties of Gaussian random fields; therefore, our results are independent 
of the shape of the power spectrum. In Sect. 2, we derive expressions for the mass, the 
angular momentum, the energy, and, consequently, for the spin parameter of a dark halo. 
In Sect. 3, the expression for the spin parameter is statistically analyzed, assuming that 
the primordial density fluctuations can be described by a Gaussian random field. The 
resulting probability distribution is then discussed in terms of the shape of the dark halo, 
the peak height, and the background cosmology. The results are then compared to the 
outcome of iV-body simulations. 



2 Fundamental parameters of dark haloes 

In this section, we derive expressions for mass, angular momentum, energy, and spin 
parameter of dark haloes within the framework of linear theory. The dimensionless spin 
parameter is defined by 



GM 5 /< 

Note that its dependence on the total energy of the system is rather weak 



2.1 Linear theory and Zel'dovich approximation 

For sufficiently early times, the evolution of density perturbations S(r,t) = [g(r,t) — 
Qb]/Qb is well described by linear theory. g(r,t) is the density at a given point (r,t) in 
spacetime, and g\j is the mean background density. The density contrast 8 fulfils Poisson's 
equation, A<£> = —8, <& being the peculiar gravitational potential. In Eulerian formulation, 
the comoving position of a perturbation in space is fixed and only its amplitude changes, 



8{t) = b(t)8(U) , 



(2) 
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where b(ti) is normalized to unity. In an Einstein-de Sitter universe, b(t) = a(t), a(t) 
being the dimensionless cosmological expansion function, normalized to unity at t = t{. 
For an open universe {Qq < 1), the linear growth factor b(t) can be approximated to 
within a few percent by 

(3) 

1 + ilz 

with 

ii= 2 - 5 "° (4) 
l + 1.5^o 

(Shandarin et al. 1983), z = (1 + z\)/a(t) — 1 being the redshift corresponding to time t. 

The accuracy of linear theory can be improved if particle trajectories are considered 
instead of the density contrast. This technique is known as Zel'dovich approximation 
(Zel'dovich 1970) and is most easily understood in the Lagrangian formulation of per- 
turbation theory (Buchert 1992). Let rj, j £ {1,2,3}, be physical particle coordinates 
starting at Lagrangian coordinates qj. Then, 

r j {q j ,t) = a{t)[q j + {b{t)-l)* tj (q)] , (5) 

where the comma followed by an index j indicates the partial derivative with respect to 
Qj- 

2.2 Angular momentum 

Consider now a fiducial trajectory rj(qj,t) which, for instance, could be the trajectory 
of the center-of-mass of the mass distribution to be considered. For a neighbouring 
trajectory, we then have 

dv ■ 

5r i = -Z- Ls Qk = a(t) [<% + (b(t) - l)<P, jk Sq k ] , (6) 
where the sum convention is implied. Since the particle velocity is 

Vj = r- = a(t) qj + j t [a(t)(b(t) - l)]$ d , (7) 
the velocity of the particle considered relative to the fiducial particle is 

S Vj = a(t)6 qj + j t [a(t)(b(t) - l)]$ jk 6q k . (8) 

Hence, the angular momentum of a volume element dV with density g at distance 5rj 
relative to the fiducial particle is 



dLj = gdVe jk i 



a(t)j t [a(t)(b(t) - l)]$ M 8q k 8q n + a(t)(b(t) - l)a(t)^ tkm 5qi5q r , 



, (9) 



where tj k \ is the totally antisymmetric tensor with €123 = 1. Integrating over the volume 
of interest T, we obtain after changing the summation indices and accounting for the 
antisymmetry of 6j k i 

Lj = a 2 (t)b(t)e jkl I kn <Pj n , (10) 
where we have defined the inertial tensor of Y 
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Ijk = dV g(ti)SqjSq k . (11) 
Jr 

This result was first derived by White (1984), who also discusses Eq.(lO) in detail. Note 
that the origin of the angular momentum is the shear flow of the cosmic medium. In 
the Zel'dovich approximation, motion is restricted to irrotational flows. This guarantees 
the existence of the potential <5 whose gradient is the velocity field and whose source 
can, in linear approximation, be identified with the density contrast field 8. Note also 
that we use the Zel'dovich approximation only for the derivation of Eq.(lO); it enters the 
derivation of the spin parameter only in this respect. 

Now, we approximate the shape of T at initial time ti by a triaxial homoeoid of 
density g with the same inertial tensor as T. Let the principal axes Uj be ordered 
such that ui > U2 > U3. We introduce the ellipticity parameter e = (^2/^1) arid the 
prolateness parameter p = (u^/ui); then, 1 > e > p. Rigorously, the Uj depend on time; 
however, since we need the inertial tensor only at ti, we omit the time dependence and 
use the abbreviation Uj = Uj(t{). After diagonalization, the inertial tensor reads 

I jk = diag(/x j ) , j e {1,2,3} , (12) 

where 

Hi = —uj , 1x2 = Hie 2 , fi 3 = ix x p 2 . (13) 
5 

The total mass M enclosed by the ellipsoid is 

M = —u\ Q ep . (14) 

Referring the density to the background density £>b of the underlying Friedmann-Lemaitre 
model, we have at an initial time t\ 

g(t i ) = g h (t i )(l + 5 i ), (15) 

where r5; is the initial density contrast of the volume. In terms of the initial cosmological 
density parameter j?;, 

«.(*.) = , (16) 

where iJ; is the Hubble parameter of the background model at t\. Combining Eqs.(14), 
(15) and (16), we obtain 

M=^(l + fc)«?cp. (17) 



From Eq.(13), we get 



and, together with Eq.(17), this yields 



= ^S*g±^!'V (19) 
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2.3 Potential energy 

A homoeoid with principal axes Uj and density g has the potential energy 

8tt 2 



W = -~—Gq 2 u\ 
15 



2eV E(M) 

sin 9 



(20) 



see, e.g., Binney & Tremaine (1987, table 2.1). In Eq.(20), F(9,k) is the incomplete 
elliptic integral of first kind, with 



e 2 

6 = arccos(p) , k = J - — — . (21) 

Together with Eq.(14), we obtain from Eq.(20) 

W = _A GM 2!M . (22 ) 
5tti smfc> 

The potential energy equals the total energy when the ellipsoid achieves its maximum 
extension. This neglects the rotational energy of the ellipsoid, which is, however, small 
compared to the potential energy, as will be confirmed below. A homogeneous sphere 
with density contrast r5; expands until it reaches the radius 



„ max „, ' u i c i /oqA 



see, e.g., Bartelmann et al. 1993, Eq.(23); there, it was also shown that the collapse 
history of a homoeoid is, at least statistically, well approximated by the collapse timescale 
of a sphere, e; = 1 — Q\ is the deviation of the initial density parameter from unity. 
That means that density fluctuations with 5\ < (3/5)e; never collapse. Therefore, it is 
convenient to introduce the quantity 

7i = - e ; , (24) 

i.e., density fluctuations with 7; > collapse, those with 7; < do not; since r5; <C 1 and 
e; 1, 7; is also much smaller than unity. Approximating the maximum value of the 
largest principal axis of the ellipsoid by (23), we obtain for the total energy 

£ = _3G 7 ,M3F(M) ( 
5iti sin 9 

Ml/2 = ^ 

Ul 

E _ 3G 7 iM 5 / 2 F(g,fe) ^ 
5 \/5a*i s i n ^ ' 

an expression for the energy which turns out convenient later on. Note that the approx- 
imation for the energy is the probable crudest in the investigation presented here. The 
resulting error, however, is relatively small because only \[E~ enters the definition of the 
spin parameter (see below). 



Moreover, from Eq.(13), 



and therefore 
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2.4 The spin parameter 

From Eqs.(lO) and (12), we obtain for the angular momentum L in the principal-axis 
system of the inertial tensor 



(28) 



j>k 

Therefore, the square of the spin parameter (1) reads 



A 
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a 2 (t)b(t) 



2 7i F(0, k) Sj>fc [(Vj ~ fik) & Jk 



H? sin 9 



(29) 



In somewhat more detail, the term F(9,k)/sm0 reflects the change in gravitational 
binding energy compared to a homogeneous sphere, i.e., it is an intrinsic property of Y. 
The rightmost part of Eq.(29) describes the tidal interaction of Y with the surrounding 
matter. Note that in the numerator the same power of [ij appears as in the denominator; 
i.e., the term depends only on the axis ratios of Y. The term <S> ^ is dimensionless and 
of the same order of magnitude as 8j. The term (a 2 b) 2 /H 2 measures the time evolution 
of the angular momentum in units of the initial Hubble parameter H{. This term will be 
considered in the following subsection. 



2.5 Dependence on time 

The time-dependent factor in square brackets in Eq.(29) remains to be evaluated. To do 
so, we start from the definition of the 'velocity factor' /, 



a b 



(30) 



see, e.g., Peebles [1993, Eqs.(5.116) and (5.120)]. With Eq.(30), the time-dependent 
factor in Eq.(29) can be written 

a 2 b = a 2 Hbf . (31) 



For sufficiently large collapse redshifts, we can employ the Einstein-de Sitter limit, 



a(t) 



2/3 



, b(t)~a(t) , U 



3#iVA 



and 



With 



it follows 



H(t) ~ Hiy/a, 



J?i = 1 - e; ~ 1 



1 - i?n 



fto(l + *i) 
a 2 (t)b(t) ~ ^HfQit ~ *H?t 9i (t) 



^(0 = ^°- 6 (t)(i- ei )^^ u - 6 W 



0.6/ 



l-^o 
f2 (l + Zi) 



(32) 

(33) 
(34) 

(35a) 
(356) 
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Fig. 1. The exponent 0(70) introduced in 
Eq.(38, solid line). The dotted line displays 
the approximation a. ~ (0. 23/70), which is 
seen to provide an excellent fit 



If a perturbation collapses early, the density parameter Q(t) at collapse time will be close 
to unity and depend only very weakly on t. Then, Eq.(35a) shows that in this case the 
spin parameter grows proportional to t as long as the linear approximation provides a 
viable description. 

We now insert for t the time when a homogeneous sphere of density contrast <5; 
achieves its maximum extension; this is given by 

7T 

-^i^max = TJo ) (36) 

2(§*- £i ) 3/2 

see, e.g., Bartelmann et al. [1993, Eq.(26)]. Inserting into Eq.(35a), we end up with 

s ?* (tm) -»/» . (s7) 

Mi 4 

It will later on turn out favourable to 'rephrase' the £ max -dependency of the function g x . 
Of course, the most straightforward procedure would be to just insert the (analytically 
known) functions a(t) and b(t) into Eq.(29) and avoid any approximations. However, 
since we are aiming on the formation of galaxies or rich clusters, we suspect the col- 
lapse time £ max to be early enough for the Einstein-de Sitter limit to be a reasonable 
approximation. Then, the further treatment can be considerably simplified. 
We start with the ansatz 

<?i(*max) = Q% (38) 

and determine the coefficient a by performing a least-square fit of the approximation (37) 
to the exact function a 2 (t)b(t). The so-determined a will depend on the amplitude of the 
density perturbation. For large overdensities, we expect a — > 0, since then the density 
parameter at the turnaround time will be very close to unity, while for low overdensities, 
a should approach 0.6. 

Figure 1 displays a as a function of 70 = 7i6(^)| z =o- The figure clearly shows that, 
in the range considered, a is very well approximated by a power law, 



(39) 



a 
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Fig. 2. The factor a 2 (t max )b(t max ) occur- 
ring in the expression (29) for the spin 
parameter, scaled by H\. The solid line 
gives the exact (numerically determined) 
value, the dotted line the approximation of 
Eq.(37) for 70 £ {1, 3, 5}, which can be con- 
sidered typical for fluctuations evolving into 
clusters of galaxies or galaxies (see text) 



This confirms that for high 70 and thus large density perturbations, a — > 0; e.g., for 
70 = 5 we have a ~ 0.05. Thus, we can write with an error of about 0.1% 

^0,70) ^o°- 23/7 ° • (40) 

Moreover, we compare in Fig. 2 the exact function a 2 (t)b(t) with the approximation 
(37), where g\ is taken from (40); curves for 70 £ {1,3,5} are displayed. The figure 
shows that the approximation (37) to the exact behaviour of the time-dependent term is 
excellent. 

Thus, the i?o dependence is rather weak, a = 0.05 ... 0.1 being typical. That means, 
that for galaxy-sized perturbations, the spin parameter should practically be independent 
of Qq. Even for much larger fluctuations (e.g., fluctuations on cluster scale), with 70 ~ 1, 
the dependence oc i?o' 23 is still weak. For consistency, we now consider g x as a function 
of e\ and r5;, g\(e\,5\) = g\ [J?n(ei), 7o(5;, e;)]. Introducing this expression into Eq.(37), the 
square of the spin parameter becomes 

2 _ 277T 2 gUeu 5Q F(9, k) E J>fc Uj - »k) 2 $% 
1000 7? sin0 Vi^WnH ' 

Note that A 2 depends on a sum of terms oc ^ 2 j k /li- In case of an Einstein-de Sitter 
universe, this term converges to a finite value for r5; — > because of Poisson's equation. 
If, however, e x > 0, it diverges for r5; — >■ |e;. This can be understood as follows: In 
an open universe, density perturbations with |r5; < e x do never collapse but 'freeze in' 
at z ~ l/^o — 1- Nevertheless, the tidal field is present, i.e., the tidal torques are 
still acting, the angular momentum grows with time. Within linear theory, such a non- 
collapsing 'dark halo' would acquire an infinite amount of angular momentum. Therefore, 
we expect that for lately formed objects, the spin parameter should be higher. On the 
other hand, the factor g\ = ^7q' 23 ^ 70 suppresses the angular momentum due to the larger 
'Hubble friction' in an open universe. This becomes more important with decreasing 
70, because for high 70, most of the gravitational evolution is finished when the cosmic 
evolution is still indistinguishable from that of a flat universe. For lower 70, a noticible 
fraction of the dynamical time is spent during the phase when Q significantly deviates 
from unity, where the 'Hubble friction' becomes more important. We will see later on 
that these two effects balance very well for objects with 70 ^ 1 (see Fig. 5). 
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Fig. 3. Contour plot of the shape-depen- 
dent function G(e,p) for (e,p) £ [0.5, 1] 



3 Statistics 



In the following section, the statistical properties of the spin parameter will be consid- 
ered within the framework of Gaussian random fields. A strictly homogeneous Gaussian 
random field is entirely determined by its power spectrum (8%) = P(k). Its variance cr;, 
given by 

d 3 k 



(27T): 



:P(k,U)W x (k) , 



(42) 



is a measure for the mean peak height at the time t[. Wx is the window function 
smoothing out fluctuations on scales smaller than the typical extension X of the con- 
sidered volume. With <7o we will design the variance at z = 0, extrapolated from z\ via 
Eq.(2). 

3.1 Ensemble average 

Averaging over ensembles, and assuming that the quantities determining the shape of 
the ellipsoid are statistically independent of <Pjk and 8[, we obtain 

(X z ) = — ( ) ( E j>k (N~ ^l x _ (43) 



1000 



7i 



sin^ ^iV^2^3 



in this equation, it was further assumed that, on average, {&,jk) is the same for all i ^ j, 
so that <Pjk can be replaced by <5 12- Hence, the spin parameter separates into a term 
dependent only on the (initial) geometry of the ellipsoid, 



{G 2 (e,p)) 



F ( arccos(p) , \ / } e l ) 

A — r 3(1 



ep^Jl — p 2 



[3(l + e 4 +/)-(l + e 2 +p 2 ) 2 ] , 



(44) 



and a term dependent on the cosmological density parameter (via e;) and the statistical 
properties of the potential which are incorporated in the perturbation spectrum, 



K 2 (<7i,ei) 



7f 



(45) 



Figure 3 displays contours of the shape-dependent function G(e,p). 



1U 
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The assumption that the shape of the volume T is statistically independent of the 
surrounding tidal field <Pjk can be motivated by the following argument. Let the tidal 
field <Pjk be decomposed into two parts <P ll j k and ^ e jj*. Sources for <P ll j k and <P e ? k are the 
matter distributions inside and outside T, respectively. The principal-axis system of <P 1Y j k 
has to be aligned with Ijk because otherwise the homoeoid would be spun up by its own 
tidal field. The alignment between the external and the internal tidal fields, if present, 
should be weak for the assumption of statistical indepence between <P e * k and $ ll ^ t k to 
be valid. iV-body simulations by Barnes & Efstathiou (1987) support this view, in that 
they exhibit only a weak tendency for the principal axes (or spin vectors) of neighbouring 
haloes to be aligned. In any case, the assumption of statistical independence tends to 
overestimate the growth of angular momentum, because if the principal-axis system of 
T tends to be partially aligned with the external tidal field, the angular momentum gain 
of T is reduced. 



3.2 Eigenvalue distribution of <&,jk 

We now statistically evaluate the cosmological factor (gf^^/li)- This evaluation is 
easiest performed rotating first into the principal-axis system of the tidal tensor <Pjk and 
applying a random rotation to the result afterwards. Averaging over the three Euler 
angles of this rotation, we will obtain the desired relation in Eq.(51) below. Since 7; 
depends only on the trace of <&jk, it is unaffected by rotations; therefore, the principal- 
axis system of gf^ 2 j k /li is identical with that of <&jk- 

Since <P^k is symmetric, there exists a local coordinate transformation T(a, (3, 7) 
dependent on the three Euler angles (0,7) G [0, 2n), (5 e [0,7r) such that 

T"V, P, 7) ($jk) T(a, (3, 7) = diagfai, <p 2 , <Ps) ■ (46) 

For a Gaussian random field <5, the probability distribution of the eigenvalues cpj of <Pjk 
is known to be 

15 3 

p((f)l,(j)2,(f)3)d 3 (f) = -= |(03 - 0i)(0 3 - 4>2){4>2 ~ 0l)| 

x exp I -- [2{<j>\ + ^1 + 4,1) - (0!0 2 + 0!0 3 + 020s)] I d 3 (f) ; 

(Doroshkevich 1970, Bartelmann & Schneider 1992). We have introduced <f>j = {(fj/o\) 
with o\ defined by the perturbation spectrum of the density contrast P(k) [Eq.(42)]. The 
initial density contrast is given by Poisson's equation, 

3 

5i = A$ = tr ($ >ifc ) = a ; ^ ^ . (48) 

i=l 

As motivated in section 3.1, the principal-axis system of the tidal tensor <Pjk can be 
assumed to be independent of the principal-axis system of the inertial tensor Ijk- Then, 
the Euler angles (a, ft, 7) between the two systems assume random values in their intervals 
of definition. Consider an arbitrary diagonal matrix Mjk = diag(mi, m-2, m^) which is 
arbitrarily rotated by a rotation matrix T(a, (3, 7), 



M' jk (m u m 2 ,m 3 ,a,/3,^) = T(a,P,j)(M jk )T 1 {a,f3, 1 ) 



(49) 
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It can then be shown by straightforward, but rather tedious algebra that, averaging over 
the Euler angles, 



(M{i)2 ) = i5 3 EK 2 )+ 2 E( m ^) > 

V i=i 3+k J 

(M[ 2 f) = 1 f ]T(m 2 ) - ^K-m,} . 

\i=i i#fc / 



(50) 



If the mj have the same statistical properties for all j, and if (rrijirik) = for j / A;, as 
is the case for the Eq.(50) can be simplified to 



( M 'n) 2 ) = f <«»?> , 
(M' jk ) 2 ) = \(m\) ; 



(51) 



the same result could have been obtained also by considering invariants of M^. For 
example, we obtain 

{4>D = 1^ <<£ 2 n> = I , fe) = ^ ; (52) 

these results agree with those given in Bardeen et al. [1986, Eq.(Al)]. 

3.3 Conditional probability 

With Eq.(51), we can write 

</i 2 (^i)f 2 12 



7; 
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/ 



d 3 (f>p{(f>i,(j)2,fo) 



(E0j - wiY 



(53) 



where u)\ = (3/5)(ei/o"i). Within this notation, the effective density contrast is given by 

The integral in Eq.(53) has to be evaluated for sufficiently large values of i.e., 
restricting to such density perturbations which collapse at moderate to high redshifts. 
We therefore choose 



(54) 



i=l 



where v ;> 1 appears as a reasonable choice. Note that because of Poisson's equation, 
A<£> = —8, the restriction (54) on the density contrast also restricts the tidal field. The 
integral in Eq.(53) then has to be replaced by 



/d 3 0p(0i,0 2 ,03) 



9, 2 (^i)^ 



H (E <t>j - Wi - I/) 



Jd 3 0p(01,^>2,03)H(5]^ -Wi-Z^) 



K 2 (ai,ei,zy) , 



(55) 



where H(x) is the Heaviside step function; Eq.(55) expresses Bayes' formula for condi- 
tional probabilities. Combining Eqs.(43), (53), and (55), we finally obtain 
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Fig. 4. The function .ff (<7;, e;, i/) for the 
Einstein-de Sitter case, e; = 0. In this 
case, K does only depend on v. The solid 
line gives K(y), the dotted line shows the 
power-law fit of Eq.(58) 



with 



{\ I ) = KK I (a u e u u) G 2 (e } p) } 
3 5 tt 2 



\2 
A 



50 3 



0.02 



(56) 
(57) 



As an example, Fig. 4 shows the function K(a[, e;, v) for the Einstein-de Sitter case, 
i.e., for e\ = 0. Then, K does no longer depend on a;, and is merely a function of v. As 
seen in the figure, K{y) is reasonably fit by a power law in the range of v given; we find 
with an error of <; 5% 



K[y) 



2.07 v - 29 



for v e [0.5... 3] . 



(58) 



For v = 1, we read off K ~ 0.48. Since G(e,p) ~ 1 for reasonable values of e and p, 
we have 

V 7 ^) ~ 0.07 , (59) 

which is in excellent agreement with numerically determined spin parameters. The small 
value of A supports the neglection of the rotational energy in Sect. 2. 3. 

3.4 Dependence on the density parameter 

We have defined e; = 1 — i?;, hence, 



fto(l + *i) ' 
cf. Eq.(34). Since a; and o"o are related by 

i + hz- x 

1 + Qz 

with £2 from Eq.(4), the quantity u>\ is given by 



(60) 



(61) 



bummary 



o 
G 

o 

b 



1 

0.8 
0.6 
0.4 

0.2 



+ I 

- 


i i I 


— 

-J 




.0. 


= 3 - 


- 




= 5 - 


- V = 1 












- v = 2 




— =- 




ill 





0.2 



0.4 



0.6 



Fig. 5. The function ^'(cto, J7o,f) for v = 
1,2 as a function of J?o. The solid (dotted, 
dashed) line is for ctq = 1 (3,5, respectively) 



where we have used z\ ^> 1. In terms of <7o and ,J2o, we write 

(a , i? , v)=K [cr; (cr , i? ) , e; (^o) , H • (63) 

In Fig. 5, we show the dependence of K' on o"o and i?o for v = 1. As expected, the 
dependency is very weak. Note, however, that for small Q$ and ctq, K' is systematically 
higher. This reflects the fact that objects collapsing late in a low-density universe are 
exposed to the tidal field for a longer period and can, therefore, acquire higher angular 
momentum. 



3.5 Probability distribution of the spin parameter 

The probability distribution for the spin parameter can also be derived in a straightfor- 
ward way. Consider the function K'(ao, 47 ,z^) introduced in Eqs.(55) and (63) as an 
average over the relevant ranges of the <f>j. Since, as demonstrated in Figs. 4 and 5, the 
dependence of K' on both v and 4?o is weak, we restrict ourselves to the Einstein-de 
Sitter case {£2$ = 1) and to v = 1. 

Then, the statistical properties of A are uniquely determined by the probability 
distribution p(4>j) for the eigenvalues of @,jk, and thus the probability distribution p(X) 
for the spin parameter can be derived using Eq.(47) and transforming variables; it is 
displayed in Fig. 6 together with the cumulative distribution 

("CO 

P(A) = dX'p(X') . (64) 

In the right panel of Fig. 6, we indicate the region where 

10% < P(A) < 90% ; (65) 

in terms of A, this interval is 

0.02 < A < 0.11 . (66) 

This is in surprisingly good agreement with results from numerical simulations. Further- 
more, the shape of the distribution function (Fig. 6, left panel) is quite similar to that 
obtained by the numerical simulations of Warren et al. (1992, Fig. 18): The probability 
distribution steeply rises for A <J 0.025, where the probability distribution reaches its 
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Fig. 6. Probability distribution (left) and 
cumulative probability distribution (right) 
for the spin parameter A in a Einstein-de 
Sitter universe with a peak height v > 1. 
The shadowed region marks the 80% inter- 
val, i.e., only 10% have lower and higher 
spin parameters, respectively. 

maximum. Then the probability begins to decrease, relatively slowly for A <; 0.1, and 
rather steeply beyond. However, the maximum and the turnaround points in the numer- 
ical simulations tend to be at lower A. Therefore, our approach tends to overestimate 
the spin parameter. Possible explanations are the neglection of the correlations between 
shape and tidal field or systematical effects resulting from our choice of v and the ellip- 
ticities (e,p), and the approximation of the turnaround time by the time of maximum 
expansion of a corresponding sphere. 

4 Summary 

We have investigated the spin parameter of dark haloes forming by gravitational collapse 
within the framework of first-order growth of angular momentum. Several approxima- 
tions enter this derivation, these are: 

1. The volume T is approximated by a homoeoid which is supposed to have the same 
inertial tensor as Y. 

2. The energy of the volume (potential and kinetic) is approximated by taking the 
potential energy of a homoeoid at maximum expansion with the axis ratio of the 
homoeoid at initial time, while this axis ratio is generally depending on time. 

3. The extension at maximum expansion of the homoeoid is approximated by the max- 
imum radius of a sphere with the density contrast of the homoeoid. 

4. The time-dependence of the spin parameter is evaluated in the high-redshift limit of 
the dynamics of the background cosmological model. 

5. The principal-axis systems of T and the tidal field are taken to be statistically 
independent. 

Assumption (5) is not obvious but seems to be well supported by numerical results. This 
assumption causes a slight overestimate of the spin parameter in the formalism above. 
This might be the reason for the slightly smaller spin parameters found in numerical 
simulations. Approximation (4) was shown in Sect. 2 to be very well satisfied. Approx- 
imations (2) and (3) only affect the energy of the homoeoid, and the spin parameter 
depends only weakly on the energy. Additionally, Bartelmann et al. (1993) showed that 
the collapse history of a homoeoid is well approximated by the collapse timescale of a 
sphere, at least in a statistical sense. Approximation (1) is not crucial because the gain 
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of angular momentum is controlled only by the inertial tensor of the mass distribution 
considered. 

On the basis of these approximations, the spin parameter can be written 

JW)- ^77^G(e,p)~0.07^-°- 29 G(e,p) , , , 

v N ' 2.07 v - 29 K ' v ' (67) 

A G [0.02 . . . 0.11] z> -0 ' 29 G(e,p) for v G [0.5, 3] . 

Its median value is A ~ 0.06, and the limits on A are its 10 and 90 percentile values. In 
Eq.(67), G(e,p) ~ 1 is determined by the shape of the dark halo, and v is its initial density 
contrast 5{ in units of a;. The agreement with numerically determined spin parameters is 
astonishing, in particular, also the probability distribution of A (Fig. 6, left panel) is well 
reproduced. The spin parameter increases linearly with time and approaches the value 
given by Eq.(67) after turnaround time, when the volume under consideration shrinks 
and makes it increasingly difficult for the external tidal field to spin it up. 

As has been discussed above, the dependence of A on the cosmological density pa- 
rameter is very weak (see Fig. 5). Although the velocity factor / ~ i? 0,6 (£) entered the 
derivation of A, most of the 'interesting' density perturbations collapse sufficiently early 
for the Einstein-de Sitter limit to be a viable approximation; then, i? is close to unity 
at turnaround time. When i? is low, i?o J> 0.2, the normalization of the density per- 
turbation spectrum today implies that perturbations which are now collapsed must have 
collapsed at higher redshift than in a high-density universe, because otherwise they could 
not have formed (see, e.g., Richstone et al. 1992, Bartelmann et al. 1993). Thus, at the 
formation epoch, i? is still close to unity even if it is low today. 

Additionally, Fig. 5 shows that the dependence of A on the variance ao of the linear 
density perturbation spectrum P(k) is weak. This is due to the fact that both the 
tidal field and the density contrast of the volume considered are linear in ao, and since 
the ratio of these quantities enters the spin parameter, it does hardly depend on the 
spectral normalization. In turn, this implies that the spin parameter of collapsing objects 
does only very weakly depend on the type of perturbation spectrum inserted, since all 
perturbation spectra are normalized by observations on the spatial scales considered here. 
Thus, A is expected to be weakly, if at all, dependent on the cosmogonic scenario; this is 
a direct consequence of our smoothing the perturbation spectrum over the mass scale of 
interest. However, additional angular momentum can be gained after turnaround time 
by accretion of lower-mass objects, whose number density, and therefore the accretion 
rate, does depend on the power spectrum. 

There is a mild dependence of A on the peak height of the perturbations considered; 
this peak height is measured in units of the spectral variance ao by the parameter v. For 
higher z>, the spin parameter decreases oc v~ - 29 . This behaviour reflects the fact that 
higher density perturbations reach their turnaround time earlier, so that their spin-up 
time is shorter. Since A is increasing linearly with time [cf. Eq.(35a)], a shorter spin- 
up time means a smaller spin parameter. This effect is mostly balanced because the 
tidal field in the vicinity of high density peaks is on average also enhanced due to the 
autocorrelation of high density peaks. Formally, this enters our derivation with condition 
(54), since, from Poisson's equation, the tidal field is restricted to high values together 
with the density contrast. Therefore, only combinations of high density fields with high 
tidal fields enter the ensemble average in Eq.(55). 
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While the numerical values for the spin parameter were previously known from 
numerical experiments, the virtue of the present derivation is that it explains the distri- 
bution of A in terms of the statistics of Gaussian random fields. Additionally, it explains 
why the spin parameter is insensitive to both i? and the underlying cosmogonic model, 
be it CDM or HDM or some cocktail derived from these. We expect that the statistics of 
the spin parameter would considerably deviate from the results given here if the initial 
density field, and thus the tidal field, would significantly deviate from Gaussian statis- 
tics. Additionally, the derivation above shows that for the spin parameters of galaxies 
or clusters of galaxies, the angular momentum gain during the linear collapse phase of 
these objects is sufficient. In a universe with non- vanishing cosmological constant, the 
only difference would be that the expansion time-scale of the universe is enlarged, which 
means that the spin-up period would be longer; then, we would expect the spin parameter 
to increase. 
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